Multi-omics resources for the Australian southern stuttering frog (Mixophyes australis) reveal assorted antimicrobial peptides.

The number of genome-level resources for non-model species continues to rapidly expand. However, frog species remain underrepresented, with up to 90% of frog genera having no genomic or transcriptomic data. Here, we assemble the first genomic and transcriptomic resources for the recently described southern stuttering frog (Mixophyes australis). The southern stuttering frog is ground-dwelling, inhabiting naturally vegetated riverbanks in south-eastern Australia. Using PacBio HiFi long-read sequencing and Hi-C scaffolding, we generated a high-quality genome assembly, with a scaffold N50 of 369.3 Mb and 95.1% of the genome contained in twelve scaffolds. Using this assembly, we identified the mitochondrial genome, and assembled six tissue-specific transcriptomes. We also bioinformatically characterised novel sequences of two families of antimicrobial peptides (AMPs) in the southern stuttering frog, the cathelicidins and ß-defensins. While traditional peptidomic approaches to peptide discovery have typically identified one or two AMPs in a frog species from skin secretions, our bioinformatic approach discovered 12 cathelicidins and two ß-defensins that were expressed in a range of tissues. We investigated the novelty of the peptides and found diverse predicted activities. Our bioinformatic approach highlights the benefits of multi-omics resources in peptide discovery and contributes valuable genomic resources in an under-represented taxon.

A chromosome-level genome assembly for the dugong (Dugong dugon).

The dugong (Dugong dugon) is a marine mammal widely distributed throughout the Indo-Pacific and the Red Sea, with a Vulnerable conservation status, and little is known about many of the more peripheral populations, some of which are thought to be close to extinction. We present a de novo high-quality genome assembly for the dugong from an individual belonging to the well-monitored Moreton Bay population in Queensland, Australia. Our assembly uses long-read PacBio HiFi sequencing and Omni-C data following the Vertebrate Genome Project pipeline to reach chromosome-level contiguity (24 chromosome-level scaffolds; 3.16 Gbp) and high completeness (97.9% complete BUSCOs). We observed relatively high genome-wide heterozygosity, which likely reflects historical population abundance before the last interglacial period, approximately 125,000 yr ago. Demographic inference suggests that dugong populations began declining as sea levels fell after the last interglacial period, likely a result of population fragmentation and habitat loss due to the exposure of seagrass meadows. We find no evidence for ongoing recent inbreeding in this individual. However, runs of homozygosity indicate some past inbreeding. Our draft genome assembly will enable range-wide assessments of genetic diversity and adaptation, facilitate effective management of dugong populations, and allow comparative genomics analyses including with other sirenians, the oldest marine mammal lineage.

Translating genomic advances into biodiversity conservation.

A key action of the new Global Biodiversity Framework is the maintenance of genetic diversity in all species to safeguard their adaptive potential. To achieve this goal, a translational mindset, which aims to convert results of basic research into direct practical benefits, needs to be applied to biodiversity conservation. Despite much discussion on the value of genomics to conservation, a disconnect between those generating genomic resources and those applying it to biodiversity management remains. As global efforts to generate reference genomes for non-model species increase, investment into practical biodiversity applications is critically important. Applications such as understanding population and multispecies diversity and longitudinal monitoring need support alongside education for policymakers on integrating the data into evidence-based decisions. Without such investment, the opportunity to revolutionize global biodiversity conservation using genomics will not be fully realized.

The genome sequence of the critically endangered Kroombit tinkerfrog ( Taudactylus pleione).

The Kroombit tinkerfrog ( Taudactylus pleione) is a stream-dwelling amphibian of the Myobatrachidae family. It is listed as Critically Endangered and is at high risk of extinction due to chytridiomycosis. Here, we provide the first genome assembly of the evolutionarily distinct Taudactylus genus. We sequenced PacBio HiFi reads to assemble a high-quality long-read genome and identified the mitochondrial genome. We also generated a global transcriptome from a tadpole to improve gene annotation. The genome was 5.52 Gb in length and consisted of 4,196 contigs with a contig N50 of 8.853 Mb and an L50 of 153. This study provides the first genomic resources for the Kroombit tinkerfrog to assist in future phylogenetic, environmental DNA, conservation breeding, and disease susceptibility studies.

Health Assessments of Koalas after Wildfire: A Temporal Comparison of Rehabilitated and Non-Rescued Resident Individuals.

Many koalas (Phascolarctos cinereus) required rehabilitation after the 2019/20 Australian megafires. Little is known about how the post-release health of rehabilitated koalas compares to non-rescued resident koalas. We evaluated health parameters in rehabilitated koalas and resident koalas in burnt and unburnt habitat in southern New South Wales, Australia. Health checks were undertaken within six weeks of fire (rehabilitated group), 5-9 months post-fire and 12-16 months post-fire. Body condition improved significantly over time in rehabilitated koalas, with similar condition between all groups at 12-16 months. Rehabilitated koalas therefore gained body condition at similar rates to koalas who remained and survived in the wild. The prevalence of Chlamydia pecorum was also similar between groups and timepoints, suggesting wildfire and rehabilitation did not exacerbate disease in this population. While there was some variation in measured serum biochemistry and haematology parameters between groups and timepoints, most were within normal reference ranges. Our findings show that koalas were generally healthy at the time of release and when recaptured nine months later. Landscapes in the Monaro region exhibiting a mosaic of burn severity can support koalas, and rehabilitated koala health is not compromised by returning them to burnt habitats 4-6 months post-fire.

Tasmanian devil cathelicidins exhibit anticancer activity against Devil Facial Tumour Disease (DFTD) cells.

The Tasmanian devil (Sarcophilus harrisii) is endangered due to the spread of Devil Facial Tumour Disease (DFTD), a contagious cancer with no current treatment options. Here we test whether seven recently characterized Tasmanian devil cathelicidins are involved in cancer regulation. We measured DFTD cell viability in vitro following incubation with each of the seven peptides and describe the effect of each on gene expression in treated cells. Four cathelicidins (Saha-CATH3, 4, 5 and 6) were toxic to DFTD cells and caused general signs of cellular stress. The most toxic peptide (Saha-CATH5) also suppressed the ERBB and YAP1/TAZ signaling pathways, both of which have been identified as important drivers of cancer proliferation. Three cathelicidins induced inflammatory pathways in DFTD cells that may potentially recruit immune cells in vivo. This study suggests that devil cathelicidins have some anti-cancer and inflammatory functions and should be explored further to determine whether they have potential as treatment leads.

Genomes of two Extinct-in-the-Wild reptiles from Christmas Island reveal distinct evolutionary histories and conservation insights.

Genomics can play important roles in biodiversity conservation, especially for Extinct-in-the-Wild species where genetic factors greatly influence risk of total extinction and probability of successful reintroductions. The Christmas Island blue-tailed skink (Cryptoblepharus egeriae) and Lister's gecko (Lepidodactylus listeri) are two endemic reptile species that went extinct in the wild shortly after the introduction of a predatory snake. After a decade of management, captive populations have expanded from 66 skinks and 43 geckos to several thousand individuals; however, little is known about patterns of genetic variation in these species. Here, we use PacBio HiFi long-read and Hi-C sequencing to generate highly contiguous reference genomes for both reptiles, including the XY chromosome pair in the skink. We then analyse patterns of genetic diversity to infer ancient demography and more recent histories of inbreeding. We observe high genome-wide heterozygosity in the skink (0.007 heterozygous sites per base-pair) and gecko (0.005), consistent with large historical population sizes. However, nearly 10% of the blue-tailed skink reference genome falls within long (>1 Mb) runs of homozygosity (ROH), resulting in homozygosity at all major histocompatibility complex (MHC) loci. In contrast, we detect a single ROH in Lister's gecko. We infer from the ROH lengths that related skinks may have established the captive populations. Despite a shared recent extinction in the wild, our results suggest important differences in these species' histories and implications for management. We show how reference genomes can contribute evolutionary and conservation insights, and we provide resources for future population-level and comparative genomic studies in reptiles.

Koala Genome Survey: An Open Data Resource to Improve Conservation Planning.

Genome sequencing is a powerful tool that can inform the management of threatened species. Koalas (Phascolarctos cinereus) are a globally recognized species that captured the hearts and minds of the world during the 2019/2020 Australian megafires. In 2022, koalas were listed as 'Endangered' in Queensland, New South Wales, and the Australian Capital Territory. Populations have declined because of various threats such as land clearing, habitat fragmentation, and disease, all of which are exacerbated by climate change. Here, we present the Koala Genome Survey, an open data resource that was developed after the Australian megafires. A systematic review conducted in 2020 demonstrated that our understanding of genomic diversity within koala populations was scant, with only a handful of SNP studies conducted. Interrogating data showed that only 6 of 49 New South Wales areas of regional koala significance had meaningful genome-wide data, with only 7 locations in Queensland with SNP data and 4 locations in Victoria. In 2021, we launched the Koala Genome Survey to generate resequenced genomes across the Australian east coast. We have publicly released 430 koala genomes (average coverage: 32.25X, range: 11.3-66.8X) on the Amazon Web Services Open Data platform to accelerate research that can inform current and future conservation planning.

Extinct in the wild: The precarious state of Earth's most threatened group of species.

Extinct in the Wild (EW) species are placed at the highest risk of extinction under the International Union for Conservation of Nature Red List, but the extent and variation in this risk have never been evaluated. Harnessing global databases of ex situ animal and plant holdings, we report on the perilous state of EW species. Most EW animal species-already compromised by their small number of founders-are maintained at population sizes far below the thresholds necessary to ensure demographic security. Most EW plant species depend on live propagation by a small number of botanic gardens, with a minority secured at seed bank institutions. Both extinctions and recoveries are possible fates for EW species. We urgently call for international effort to enable the latter.

Best genome sequencing strategies for annotation of complex immune gene families in wildlife.

BACKGROUND: The biodiversity crisis and increasing impact of wildlife disease on animal and human health provides impetus for studying immune genes in wildlife. Despite the recent boom in genomes for wildlife species, immune genes are poorly annotated in nonmodel species owing to their high level of polymorphism and complex genomic organisation. Our research over the past decade and a half on Tasmanian devils and koalas highlights the importance of genomics and accurate immune annotations to investigate disease in wildlife. Given this, we have increasingly been asked the minimum levels of genome quality required to effectively annotate immune genes in order to study immunogenetic diversity. Here we set out to answer this question by manually annotating immune genes in 5 marsupial genomes and 1 monotreme genome to determine the impact of sequencing data type, assembly quality, and automated annotation on accurate immune annotation. RESULTS: Genome quality is directly linked to our ability to annotate complex immune gene families, with long reads and scaffolding technologies required to reassemble immune gene clusters and elucidate evolution, organisation, and true gene content of the immune repertoire. Draft-quality genomes generated from short reads with HiC or 10× Chromium linked reads were unable to achieve this. Despite mammalian BUSCOv5 scores of up to 94.1% amongst the 6 genomes, automated annotation pipelines incorrectly annotated up to 59% of manually annotated immune genes regardless of assembly quality or method of automated annotation. CONCLUSIONS: Our results demonstrate that long reads and scaffolding technologies, alongside manual annotation, are required to accurately study the immune gene repertoire of wildlife species.

Restoring faith in conservation action: Maintaining wild genetic diversity through the Tasmanian devil insurance program.

Conservation breeding programs aim to maintain 90% wild genetic diversity, but rarely assess functional diversity. Here, we compare both genome-wide and functional diversity (in over 500 genes) of Tasmanian devils (Sarcophilus harrisii) within the insurance metapopulation and across the species' range (64,519 km2). Populations have declined by 80% since 1996 due to a contagious cancer, devil facial tumor disease (DFTD). However, predicted local extinctions have not occurred. Recent suggestions of selection for "resistance" alleles in the wild precipitated concerns that insurance population devils may be unsuitable for translocations. Using 830 wild samples collected at 31 locations between 2012 and 2021, and 553 insurance metapopulation devils, we show that the insurance metapopulation is representative of current wild genetic diversity. Allele frequencies at DFTD-associated loci were not substantially different between captive and wild devils. Methods presented here are valuable for others investigating evolutionary potential in threatened species, particularly ones under significant selective pressures.

A targeted approach to investigating immune genes of an iconic Australian marsupial.

Disease is a contributing factor to the decline of wildlife populations across the globe. Koalas, iconic yet declining Australian marsupials, are predominantly impacted by two pathogens, Chlamydia and koala retrovirus. Chlamydia is an obligate intracellular bacterium and one of the most widespread sexually transmitted infections in humans worldwide. In koalas, Chlamydia infections can present as asymptomatic or can cause a range of ocular and urogenital disease signs, such as conjunctivitis, cystitis and infertility. In this study, we looked at differences in response to Chlamydia in two northern populations of koalas using a targeted gene sequencing of 1209 immune genes in addition to genome-wide reduced representation data. We identified two MHC Class I genes associated with Chlamydia disease progression as well as 25 single nucleotide polymorphisms across 17 genes that were associated with resolution of Chlamydia infection. These genes are involved in the innate immune response (TLR5) and defence (TLR5, IFNγ, SERPINE1, STAT2 and STX4). This study deepens our understanding of the role that genetics plays in disease progression in koalas and leads into future work that will use whole genome resequencing of a larger sample set to investigate in greater detail regions identified in this study. Elucidation of the role of host genetics in disease progression and resolution in koalas will directly contribute to better design of Chlamydia vaccines and management of koala populations which have recently been listed as "endangered."

DNA metabarcoding reveals a broad dietary range for Tasmanian devils introduced to a naive ecosystem.

Top carnivores are essential for maintaining ecosystem stability and biodiversity. Yet, carnivores are declining globally and current in situ threat mitigations cannot halt population declines. As such, translocations of carnivores to historic sites or those outside the species' native range are becoming increasingly common. As carnivores are likely to impact herbivore and small predator populations, understanding how carnivores interact within an ecosystem following translocation is necessary to inform potential remedial management and future translocations. Dietary analyses provide a preliminary assessment of the direct influence of translocated carnivores on a recipient ecosystem. We used a metabarcoding approach to quantify the diet of Tasmanian devils introduced to Maria Island, Tasmania, a site outside the species' native range. We extracted DNA from 96 scats and used a universal primer set targeting the vertebrate 12S rRNA gene to identify diet items. Tasmanian devils on Maria Island had an eclectic diet, with 63 consumed taxa identified. Cat DNA was detected in 14% of scats, providing the first instance of cats appearing as part of Tasmanian devil diets either via predation or scavenging. Short-tail shearwaters and little penguins were commonly consumed, corresponding with previous surveys showing sharp population declines in these species since the introduction of Tasmanian devils. Our results indicate that the introduction of carnivores to novel ecosystems can be very successful for the focal species, but that commonly consumed species should be closely monitored to identify any vulnerable species in need of remedial management.

Epigenetic clock and methylation studies in marsupials: opossums, Tasmanian devils, kangaroos, and wallabies.

The opossum (Monodelphis domestica), with its sequenced genome, ease of laboratory care and experimental manipulation, and unique biology, is the most used laboratory marsupial. Using the mammalian methylation array, we generated DNA methylation data from n = 100 opossum samples from the ear, liver, and tail. We contrasted postnatal development and later aging effects in the opossum methylome with those in mouse (Mus musculus, C57BL/6 J strain) and other marsupial species such as Tasmanian devil, kangaroos, and wallabies. While the opossum methylome is similar to that of mouse during postnatal development, it is distinct from that shared by other mammals when it comes to the age-related gain of methylation at target sites of polycomb repressive complex 2. Our immunohistochemical staining results provide additional support for the hypothesis that PRC2 activity increases with later aging in mouse tissues but remains constant in opossum tissues. We present several epigenetic clocks for opossums that are distinguished by their compatibility with tissue type (pan-tissue and blood clock) and species (opossum and human). Two dual-species human-opossum pan-tissue clocks accurately measure chronological age and relative age, respectively. The human-opossum epigenetic clocks are expected to provide a significant boost to the attractiveness of opossum as a biological model. Additional epigenetic clocks for Tasmanian devil, red kangaroos and other species of the genus Macropus may aid species conservation efforts.

Incomplete lineage sorting and phenotypic evolution in marsupials.

Incomplete lineage sorting (ILS) makes ancestral genetic polymorphisms persist during rapid speciation events, inducing incongruences between gene trees and species trees. ILS has complicated phylogenetic inference in many lineages, including hominids. However, we lack empirical evidence that ILS leads to incongruent phenotypic variation. Here, we performed phylogenomic analyses to show that the South American monito del monte is the sister lineage of all Australian marsupials, although over 31% of its genome is closer to the Diprotodontia than to other Australian groups due to ILS during ancient radiation. Pervasive conflicting phylogenetic signals across the whole genome are consistent with some of the morphological variation among extant marsupials. We detected hundreds of genes that experienced stochastic fixation during ILS, encoding the same amino acids in non-sister species. Using functional experiments, we confirm how ILS may have directly contributed to hemiplasy in morphological traits that were established during rapid marsupial speciation ca. 60 mya.

Future-proofing the koala: Synergising genomic and environmental data for effective species management.

Climatic and evolutionary processes are inextricably linked to conservation. Avoiding extinction in rapidly changing environments often depends upon a species' capacity to adapt in the face of extreme selective pressures. Here, we employed exon capture and high-throughput next-generation sequencing to investigate the mechanisms underlying population structure and adaptive genetic variation in the koala (Phascolarctos cinereus), an iconic Australian marsupial that represents a unique conservation challenge because it is not uniformly threatened across its range. An examination of 250 specimens representing 91 wild source locations revealed that five major genetic clusters currently exist on a continental scale. The initial divergence of these clusters appears to have been concordant with the Mid-Brunhes Transition (~430 to 300 kya), a major climatic reorganisation that increased the amplitude of Pleistocene glacial-interglacial cycles. While signatures of polygenic selection and environmental adaptation were detected, strong evidence for repeated, climate-associated range contractions and demographic bottleneck events suggests that geographically isolated refugia may have played a more significant role in the survival of the koala through the Pleistocene glaciation than in situ adaptation. Consequently, the conservation of genome-wide genetic variation must be aligned with the protection of core koala habitat to increase the resilience of vulnerable populations to accelerating anthropogenic threats. Finally, we propose that the five major genetic clusters identified in this study should be accounted for in future koala conservation efforts (e.g., guiding translocations), as existing management divisions in the states of Queensland and New South Wales do not reflect historic or contemporary population structure.

Threatened Species Initiative: Empowering conservation action using genomic resources.

Globally, 15,521 animal species are listed as threatened by the International Union for the Conservation of Nature, and of these less than 3% have genomic resources that can inform conservation management. To combat this, global genome initiatives are developing genomic resources, yet production of a reference genome alone does not conserve a species. The reference genome allows us to develop a suite of tools to understand both genome-wide and functional diversity within and between species. Conservation practitioners can use these tools to inform their decision-making. But, at present there is an implementation gap between the release of genome information and the use of genomic data in applied conservation by conservation practitioners. In May 2020, we launched the Threatened Species Initiative and brought a consortium of genome biologists, population biologists, bioinformaticians, population geneticists, and ecologists together with conservation agencies across Australia, including government, zoos, and nongovernment organizations. Our objective is to create a foundation of genomic data to advance our understanding of key Australian threatened species, and ultimately empower conservation practitioners to access and apply genomic data to their decision-making processes through a web-based portal. Currently, we are developing genomic resources for 61 threatened species from a range of taxa, across Australia, with more than 130 collaborators from government, academia, and conservation organizations. Developed in direct consultation with government threatened-species managers and other conservation practitioners, herein we present our framework for meeting their needs and our systematic approach to integrating genomics into threatened species recovery.

Genome assembly of the numbat (Myrmecobius fasciatus), the only termitivorous marsupial.

The numbat (Myrmecobius fasciatus) is an endangered Australian marsupial, and the last surviving member of the Myrmecobiidae family. The numbat regularly undergoes torpor and is unique amongst marsupials as it is the only diurnal and termitivorous species. Here we sequenced the first draft genome of the numbat using 10× Genomics Chromium linked-read technology, resulting in a 3.42 Gbp genome with a scaffold N50 of 223 kbp. A global transcriptome from liver, lung and tongue was also generated to aid genome annotation, identifying 21,465 protein-coding genes. To investigate adaptation to the numbat's termitivorous diet and arid/semi-arid range, we interrogated the most highly expressed transcripts within the tongue and manually annotated taste, vomeronasal and aquaporin gene families. Antimicrobial proteins and proteins involved in digestion were highly expressed in the tongue, alongside umami taste receptors. However, sweet taste receptors were not expressed in this tissue, which combined with the putative contraction of the bitter taste receptor gene repertoire in the numbat genome, may indicate a potential evolutionary adaptation to their specialised termitivorous diet. Vomeronasal and aquaporin gene repertoires were similar to other marsupials. The draft numbat genome is a valuable tool for conservation and can be applied to population genetics/genomics studies and to investigate the unique biology of this interesting species.